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The transition form factor for electrodisintegration of a two-body bound system is calculated 
in the Bethe-Salpeter framework. For the initial (bound) and the final (scattering) states, we use 
our solutions of the Bethe-Salpeter equation in Minkowski space which were first obtained recently. 
The gauge invariance, which manifests itself in the conservation of the transition electromagnetic 
current J ■ q — 0, is studied numerically. It results from a cancellation between the plane wave and 
the final state interaction contributions. This cancellation takes place only if the initial bound state 
BS amplitude, the final scattering state and the operator of electromagnetic current are strictly 
consistent with each other, that is if they are found in the same dynamical framework. A reliable 
result for the transition form factor can be obtained in this case only. 
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I. INTRODUCTION 


Computing the electromagnetic (EM) form factors in the Bethe-Salpeter (BS) approach [l| requires the 
solutions of the BS equation in Minkowski space. 

The main reason is that the Wick rotation 0 , which allows to go from Minkowski to Euclidean space in the 
BS equation, cannot be performed in the integral expression of the EM form factor (see e.g. i). However, in 
contrast to the Euclidean case, finding the Minkowski space solution is complicated by the many singularities in 
the integrand of the BS equation and in the amplitude itself. In the recent years, and using different independent 
methods, these difficulties have been overcome and an important progress was achieved. 

In one of these methods [^, the kernel of the BS equation is approximately represented in a separable form. 
This allows to considerably advance analytically and therefore simplifies finding the solution. 

In the method developed in refs. the BS amplitude is represented as an integral over a weight function g 
- the so called Nakanishi transform [7| - which satisfies a nonsingular equation. A modification of this method, 
based on the light-front projection of the BS amplitude, was developed in [iSl and used to find the bound 
state Minkowski BS amplitude. The elastic EM form factor was also calculated in [l^. An equation for the 
Nakanishi function for the scattering states was derived [T^ . 



FIG. 1: Left panel: Feynman graph for contribution of FSI to the transition EM form factor. Right panel: PW 
contribution to the transition EM form factor. 


Another method [illii is based on the direct solution in Minkowski space of the BS equation after an 
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appropriate treatment of singularities. The scattering problem was there solved and the off-mass shell scattering 
amplitude first computed. Once reduced to the mass shell, the latter reproduces the phase shifts. Taken off-mass 
shell, it allows to calculate the electrodisintegration of the bound system, i.e. the form factor of the transition 
bound to scattering state. 

An important contribution to this form factor, incorporating the final state interaction (FSI), is given by the 
Feynman graph shown in Fig. [T] (left panel). The right and left vertices in this graph are just the Minkowski BS 
amplitudes for the bound (left) and scattering (right) states. If both vertices correspond to a bound state (case 
of the elastic form factor), the Nakanishi transform allows to calculate the 4D Feynman integral corresponding 
to Fig. [T](left panel), with an integrand containing three singular propagators, analytically [l3|. Then the non¬ 
singular integral with the weight functions g is safely calculated numerically. For the scattering state, though 
the Nakanishi transform also exists [l^, the corresponding weight function g at positive energies is not yet 
computed. Note, however, that very recently g was found in the zero-energy limit that allowed to calculate the 
scattering length IT^ . Without using the Nakanishi representation the scattering state vertex can be obtained 
only numerically |15l | and therefore the singular 4D Feynman integral corresponding to Fig. [T](left panel) must 
be computed numerically as well. This calculation, providing the transition electromagnetic current and the 
form factor, requires however some care to take properly into account the pole singularities of the propagators. 

The aim of this paper is to give the detail of the first results presented in and analyze the conservation 

of the calculated electromagnetic current in the inelastic transition. We will see that this current is indeed 
conserved, as it should be from general principles [2^. However this conservation is due to a rather delicate 
cancellation between the plane wave (PW) contribution (right panel of Fig. [T|) and the final state interaction 
(left panel of Fig. [T]) which requires a strict consistency between, on one hand, the bound and scattering state 
solutions and, on the other hand, the electromagnetic current operator. It thus provides a strong test for all 
these quantities simultaneously. 

The need for an internal consistency between states, currents and dynamical equation to ensure the gauge 
invariance was extensively discussed in in the framework of the BS and the Gross spectator equations. Our 
numerical results are in agreement with this general expectation. We will show that if this consistency and 
hence the gauge invariance is violated, a consequence of that is not only the appearance of a non-conserved part 
in the current - which anyway drops out in the cross section - but that the current as a whole is not valid at 
all. In other words, the transition form factors extracted from the conserved part of the non-conserved current, 
are also deficient. 

In order to illustrate our treatment of the singularities, we will restrict to the spinless particles. The gen¬ 
eralization to the fermion case is straightforward since the fermion and scalar propagators have the same 
singularities. 

The paper is organized as follows. In Section [IT] we discuss the decomposition of the transition current in the 
form factors without assuming the current conservation. In Sections IIIII and lIVi the FSI and PW contributions 
in the conserved part of the current are calculated. Section El is devoted to the discussion of the current 
conservation. In particular, the FSI and PW contributions of the non-conserved part of the current have been 
calculated and we have shown that they cancel each other. Some selected numerical results are presented in 
Section IVTl Section IVlII contains the concluding remarks. The cumbersome details of the calculations are given 
in the Appendices El [B] and [C] 


II. TRANSITION FORM FACTOR 

In the case of spinless particles, and without supposing the current conservation, the general form of the 
electromagnetic current involves two form factors: 

^ = {Pf.+P^)F,{Q^) + {p'^-p^)F2{Q^) ( 1 ) 

The decomposition dD), together with the scalar character of the constituents, implies that the initial and final 
states have total zero angular momenta, i.e. that they are composed of S-waves only. 
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In order to study the current conservation, it is convenient to redefine the form factors by introducing the 
following linear combinations F and F': 


F = Fi Fi = F 

F' = Fi-^F2 ^ F2 = §I(F-F') 

with 

9m = P^-Pm 

Ql = M'^ - (2) 


M is the initial bound state mass, M' is the invariant mass of the final scattering state. In terms of them, the 
current o can be rewritten in the form: 


J,= 


(Pm +Pm) + (Pm “Pm)^ 


W)-(p;.-pm)^W) 


( 3 ) 


Since 


q-J = QlF'{Q^) (4) 

the current conservation g • J = 0 is equivalent to F'(g^) = 0. 

Notice that in the elastic case the form factor F' is absent since the term ^ (p'^ — in ([3]) is forbidden 

by the symmetry between initial and final states. 

Notice also that the form factor F'{Q^), even if it is not zero, does not contribute to the electrodisintegration 
amplitude A. Indeed, this amplitude is given by: 

It contains the electron spinors u{k) and u{k'). Substituting here the current ([3]) and using the Dirac equation, 
we see that the term containing F'{Q'^) drops out since 

{p'^ - Pi^)Hf^'h^u{k) = u{k'){^- }i')u{k) = 0. 

Below we will calculate each of these form factors - F and F' - as a, sum of FSI (left panel in Fig. [T]) and 
PW (right panel in Fig. [T|) contributions, in the form: 

F^neliQ^) = Ff,,{Q^)+Fp^{Q^) 

FLeim = + F;^{Q^) (5) 

We will check that the full current is conserved, that is, for any the contributions to Fi„e/(0^) of 
and PW cancel each other: 


FLeiiQ^) = F'fAQ") + f ;„( q 2 ) = 0 , ( 6 ) 

provided the bound and scattering states are solutions of the BS equation with the one-boson exchange kernel. 
In this case, the EM current corresponding to the interaction of a photon with a constituent is free. We are 
however interested in a quantitative measure of the accuracy of this cancellation in a real calculation. This is 
the reason for introducing in (jS]) the non-conserved part - proportional to {p'^ — Pp)F'(Q‘^) - and the value of 
the form factor F'{Q^) will give us this measure. 

We will calculate separately the FSI and PW contributions to the form factor F'(Cp) and see with what 
accuracy they cancel each other in the sum (|6|). 
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III. FINAL STATE INTERACTION 


We start with the FSI contribution. It is obtained W applying the Feynman rules to the left panel graph of 
Fig. [Hand has the form (following the convention of [21|): 




d^k - 2fc^)r, {\p-k,p) Tf {^p' - k,p') 

(27r)4 + ie)[{p — kY — + ie][{p' — kY — + ie] ’ 


(7) 


Here F^ is the initial (bound state) vertex and F/ is the final vertex (half-off-shell scattering BS amplitude). 
As mentioned, both vertex functions were found numerically by solving the S-wave BS equation in [l5|. More 
precisely the function Fy is related to the scattering wave solution Fq by Eq. (IB29I) from Appendix IB 41 

The integrals of the type o are usually calculated by applying to the product of propagators the Feyn¬ 
man parametrization and then performing the Wick rotation. However, besides the product of propagators, 
expression © contains the initial (F^) and final (F^) BS amplitudes which are known numerically. Therefore 
the Feynman parametrization cannot be applied and we should calculate this 4D singular integral numerically, 
though after some transformations. 


It is convenient to carry out this calculations in the system of reference where p'g = po (i.e. qo = 0) and p 
and p' are collinear, i.e. they either are parallel or anti-parallel to each other, depending on the kinematical 
conditions. In the elastic case it coincides with the Breit frame p + p' = 0 and, one has of course = |p'|, 
p'o = Po- In the inelastic case, in the frame with Pq = pq we have \p\ Y \p'\- Some useful kinematical relations 
valid in this reference system are given in Appendix [^ 

In this reference frame we take the zero component of the current ([3]) and get the relation: 

Jo = 2poE(Q2) (8) 


That is: 

p _ ■ f dko (Pk (po - fco) hi {\p - k,p) Tf (ip' - k,p') 

/-WI V ( 27 r )4 PO iP^-el + ze)[ipo-koY-el_j: + ze][ipo-koY-e},_^^+ie] 

with = \/m? -I- (Y and similar expressions for and obtained using (IA4I) and (IA5I) . 


As detailed in the Appendix IB 41 in case of initial and final S-waves, all kinematical variables as well as the 
arguments of the vertex functions F appearing in ([^ can be expressed in terms of \p\, \p'\ and the integration 
variables (fcp, 2 :, |fc|) with z = k ■ p. To lighten the writing we will denote hereafter abusively p = |p|, p' = |p'| 
and k = \k\. 

After a trivial integration over the azimutal angle, the integration measure in @ becomes 

dko d^k dko dz Pdk 
(27r)'^ (271")^ 


Let us introduce the following notations, making explicit only the dependence on the integration variables: 

G{ko,z,k) 


f{ko,z,k) = 


{*0 - 4 + *e) [(Po -koY- £l_i; + *e] [{Po - koY - + *e] ’ 


p' — k 


where 


G(fco, z, k) = Ti (^^P - T/ i^^P' - k,p'^ ■ 


( 10 ) 


( 11 ) 
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Each pole singularity in (1101) is represented as a sum of its principal value and a delta-function and therefore 
the function / takes the form: 


= G{ko,z,k) 


X 

PV , 



[Po 

X 

PV 



[Po 

= /a + /2 + fi, 


1 


fcn - 


- iirSikl - e|) 


c2_ ^ 

p—k 


si 

p' — k 


- iTrS(^{po - ko f - 

- iTTS(^{po - kof - 


( 12 ) 


where /a is the contribution of the product of three principal values and no delta-functions (one single term), 
/2 is the contribution of the product of two principal values and one delta-function (three terms) and fi is the 
contribution of the product of one principal value and two delta-functions (also three terms). The product of 
three delta-functions does not contribute since their arguments cannot be zero simultaneously. 

These functions fi have the following explicit form: 


f3{ko,z,k) = G{ko,k)FY 


kn - . 


PV 


(po - koY - £ 


PV 


1 


p—k 


{po-koY - £?., 


(13) 


p' — k 


f2iko,z,k) = -inSikl-£l)G{ko,z,k)PY 


PV- 


(Po - koY - si - {po - koY - s% - 

p—k p'—k 




-ittSYpo - kof - £ 2 _,_ )G(fco, z, k) PV 


1 


PV- 


1 


’ kl-el {po-kof-s^ 


(14) 


p—k 


fi{ko,z,k) = -Tr^6{k'^-£l)S{{po-kof-e'i_j:)G{ko,z,k)PY 


{po -kof- ej,_- 


-TT S{ko - e^fSiipo - ko) - Sp_^)Giko, z, k) PYj-—^-^^^ 


p—k 


-TT^Sfpo - kof - el^Sfpo - kof - £^-,_f)G{ko, z, k) PV ^ ^ ^ (15) 

The index of fi {i = 1, 2, 3) denotes the number of the principal value products that involves. 

Our task now is to calculate the 4D integral ([9]), rewritten as 

-?^/s*(( 5 ^) = J dkodzk^dk {f3{ko,z,k) + f2{ko,z,k) + fi{ko,z,k)} 

with fi given by Eq. (11311151) . Part of this integration is calculated analytically and the remaining part, once 
transformed into a non-singular integrand, numerically. 

For calculating the singular principal value integrals in /a we will use the subtraction technique. That is, we 
subtract and add to /a an appropriately chosen singular function ho which, in variable ko has the same poles 
as /a and has no any other singularities: 


/s = (/a — h.3) + ho. 


In the difference (/a — ho), the pole singularities cancel each other and the result is a smooth function, whereas 
in the additional term ho the integral over dko is calculated analytically. 
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After this calculation, there still remains a singular expression in variable k. It is however logarithmic and 
can be treated by using standard numerical techniques, like variable change or by simply increasing the number 
of integration points. The details of all these calculations are given in Appendix [BJ 

In the integrals containing the functions /2 and /i, the integration over fcp is easily performed analytically 
by means of the delta-functions. After that, and a trivial azimuthal integration, the result is reduced to a two- 
and one-dimensional numerical integrations respectively. 

The final result for the FSI contribution © reads: 

FfsiiQ"^) = jdkodzk'^dk {/^{ko, z,k) + f 2 iko, z,k) + fi{ko, z,k)} 

= F^iQ^) + F2{Q^) + F,{Q^) (16) 

where Fi{Q'^) are defined in Appendix [Bl bv Eqs. (IB9I) . (IB13|) and (IB26|) . 


IV. PLANE WAVE CONTRIBUTION 


This contribution is displayed in the right panel in Fig. [TJ According to the Feynman rules it reads: 

P' 


d^,pW - 


[p — ky — + ie \ 2 ' 


(17) 


The delta-function follows from the four-momenta conservation in Fig. (TJ right panel: 


5^'^\k-pi) = [k-ps-^ 


We have introduced the total p' and relative Ps four-momentum of the final (non-interacting) particles 


2ps = Pl-P2 
p' = Pi + P2 


The spatial part of Ps in the rest frame of the final system p' = 0, determines the invariant final state mass 
M' = 2 w? + ps'^. One could calculate the integral over d'^k by means of the delta-function. It is however 
interesting to keep this delta-function and carry out the integration later, once extracted the S-wave from the 
final state. 

Like in the case of FSI, the form factor can be found by applying Eq. ([5]) to the Jo.pu, component of Eq. (113, 
in the system of reference where g'o = 0- That is: 


p 

± n- 


{Po - fcp) Fj (I - k,p) 

Po [{p — ky — nF ie] 


d^pB j(4) 
47r 



d'^k 


(18) 


We have introduced here the additional integration over in the rest frame p' = 0 of the final state. 

We remind that in the FSI contribution, calculated in the previous section, we decomposed the final state BS 
amplitude Fy in partial waves and took into account the S-wave only. The delta-function in (1181) replaces now 
the final BS amplitude F/. Averaging this delta-function over the solid angle Ps in the rest frame p' = 0 allows 
to select the partial S-wave in the plane wave. This is the meaning of the integral over in (flSl) . 

The integration over dQp^ and part of the integration over d'^k in (1181) are done analytically in Appendix [Cl 
The final result reads: 


Fpw — 


(po - fcp) 
Po 




1 M'k 

(jj — ky — irF -\- it 2ekPsP' 


dk 


(19) 
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where the integration limits fcip are defined in (IB15I) and Ps = \JjA, — m?. In expression (ITOl) one must 

insert /cq = \/and, in the scalar product k ■ p = kopo — zkp take the value z = zq given by Eq. (IB14I) . 
Variables po and p are the components of the initial four-momentum. The value p' is the spatial part of the 
total final state four-momentum in the frame where Pq = Po- All these components are expressed in terms of 
the momentum transfer in Eqs. (IA8I1 and (jA9|) from Appendix lAl 

Let us precise the arguments of E^ in m- They are defined analogously to the case of ESI, Eq. 

(jB27p in Appendix IB 41 Namely, solving the BS equation in the rest frame p = 0, we find ri(fco, |fc|) , where fcg 
and \k\ are also defined in the rest frame p = 0. We should express them in the frame where qo = 0. These 
expressions are given in Appendix IB 41 That is, we have to insert in (11911 the function ri(fco, |fc|) with arguments 
fco and \k\ given by the first line of Eq. (IB28I1 from the Appendix IB 41 

h = - ^{kopo- kpz), 


1^1 = \l^{koPQ-kpz)‘^-kl + P 


with fco = y/m^ + k"^ and zo defined in (IB14I1 . 


To summarize these last two sections we would like to emphasize that: (i) the full PW contribution in the 
current is given by the simple equation (HID, (a) the expression (fT^ corresponds to the S-wave projection of the 
final plane wave state, and (in) the full transition form factor - including both the ESI and PW contributions 
- is given by the sum ([S|), with Ffsi determined by Eq. (ITBl) and Fp^ by (fTOj) . 


V. CURRENT CONSERVATION 

As follows from Eq. o, the conservation of the electromagnetic current q - J = 0 implies F'{Q^) = 0 for any 
value of that is F'{Q‘^) = 0. 

To ensure this conservation, all the contributions to the current, containing the interaction of a photon with 
a charged particle, must be taken into account. In other words, in an interacting system, the true current is, in 
general, not the free one. 

For example, for the kernel given by the sum of ladder and cross-ladder, the full EM current should contain, 
in addition to the two (free) contributions displayed in Fig. [1] the cross-ladder ESI contribution shown in Fig. [5] 
and similar cross-ladder contribution for the plane wave (we suppose the exchanged particle to be neutral). The 
latter contributions are not free: they contain the interaction of constituents. The sum of four contributions 
- Fig. [T] (left and right panels). Fig. [5] (cross-ladder with ESI) and the corresponding cross-ladder PW (not 
shown) - must be conserved. 

In the case of the ladder kernel, the diagrams displayed in Fig. [T]provide the only contributions to the current. 
Therefore, the current determined by these two graphs has to be conserved if the initial and final BS amplitudes 
are also obtained with the ladder kernel. At the same time, the expressions for the contributions o and HU) 
to the current in terms of the BS amplitudes are universal - they are the same for any BS amplitude (found 
with any kernel). The conservation of their sum is provided by the particular properties of the BS amplitudes 
determined by the ladder kernel. The current is not conserved, if in Eqs. dZD and CZl) one substitutes other BS 
amplitudes (not the ladder ones). Therefore the current conservation (if any) provides a very strong test for 
the solutions themselves. In the this section we will calculate the form factor F'{Q^) and in the next section 
we will check numerically, whether it is identically zero or not. From (j4]) it follows that 

J ■ q 


F'iQ^) 


Ql 


( 20 ) 

















q 



} 


=p' 


FIG. 2: Cross-ladder contribution to the transition EM form factor. 


In the expression ([7]) for J/si, after multiplying it by consider, for a moment, only the factor 


■^{p' - p)iP + P' - 2fc) = -M^ - 2{p' -p)-k 


= 1-2 


Po=PO 


\/^zk \ 

j 


( 21 ) 


Instead of the function (^(fco, z, k) defined in Eq. (fTT|l and given by Eq. (IB27P in Appendix IB 41 we introduce 
the function: 


G'iko,z,k) 


1 - 2 


\/^zk \ 

Ql J 


ri(ko, k)rf(k'o, k'). 


( 22 ) 


G'{ko, z, k) is obtained from G{kQ,z, k), replacing the factor {po — ko)/po in G{ko,z, k) by the factor ((711) : 


(Po - fcp) 
Po 


1 - 2 


s/^zk\ 

Ql J’ 


(23) 


where 

initial 

The 


as always z is the cosine of the angle between the integration variable k and the momentum p of the 
(bound state) system in the reference frame where go = 0. 

ESI contribution to the full form factor 


F'iQ^)=F},,{Q^) + F;^iQ^) 


(24) 


is given by the same formulas than for Ffsi{Q^), eq. (ITBl) . with the replacement G{ko, z, k) —>■ G"(fco, z, k). 

The PW contribution Fp^{Q^) is calculated in a similar way. Namely Fp^{Q'^) is given by Eq. (fT^ with the 
replacement (1751) . The value fco is the same used in Eq. (ITT)) and zq is defined in Eq. (IB14I) from Appendix IB 21 
In order to obtain F'{Q‘^) identically zero the two contributions ESI and PW to the full form factor F'{Q‘^) Eq. 
(ITT)) must cancel each other. Numerically this condition is never fulfilled exactly. The value of F'(Q^) will be 
rather compared to F{Q'^) and the current conservation would manifest itself in the fact that F'{Q'^) << F{Q'^) 
for any value of Q^. 
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FIG. 3: (color online). Transition EM form factor F{Q^) as a function of . Initial (bound) state corresponds to the 
binding energy B = 0.01m; final (scattering) state corresponds to a relative momentum pa = 0.1m (final state mass 
M' = 2.00998 m). FSI contribution is shown by the dashed curve and the PW one by a dotted curve. Full form factor 
is shown by the solid curve. Left panel is the real part of form factor and right panel is the imaginary part. 


VI. NUMERICAL RESULTS 


As an example, we have calculated the transition form factor for the initial (bound) state binding energy 
B = 0.01 m (initial state mass M = 1.99 m) and for two values of the final (scattering) state relative momentum 
Ps = 0.1m and Ps = 0.5 m with corresponding final state masses M' = 2yjw? + pi values M' k. 2.00998m and 
M' « 2.236m. 



FIG. 4: (color online). The same as on the left panel of Fig. [3] for the tail of form factor 1 < < 5. 


In contrast to the elastic scattering, the inelastic transition form factor is complex. Its real and imaginary 
parts as a function of for ps = 0.1 m are shown in Fig. [31 at left and right panels correspondingly. One can 
see that at relatively small momentum transfer < 1 both contributions - FSI and PW - are important and 
they considerably cancel each other. 

The tail of the real part of form factor for > 1 and ps = 0.1m is shown in Fig. [H In this momentum 
region, FSI dominates, especially when increases. 
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FIG. 5: (color online). The same as in Fig. [3]for the final relative momentum ps = 0.5 m (final state mass: M' = 2.236 m). 


This is a natural behavior for the kinematics corresponding to Fig. IT] Indeed, due to small binding energy 
{B = 2m — M = 0.01 to) the constituents in the initial state have small relative momentum. In the scattering 
process the photon transfers the large value to one of the constituents only. However, since their relative 
energy in the final state is also small {M' — 2m ~ 0.01 to), both constituents have also small relative momentum. 
Therefore they move practically in the same direction, having both large total momentum. Since the second 
constituent does not interact with the photon, it can obtain a large total momentum only due to a strong 
interaction with the first constituent. This explains why in this kinematics the final state interaction (re¬ 
scattering) determines the tail of the form factor and dominates over the plane wave. 



q’ 



FIG. 6: (color online). Transition EM form factor F'{Q^) as a function of for ps = 0.1. Other parameters and 
notations are the same as in Fig. |3] 


The transition EM form factor for larger final relative momentum ps = 0.5 to, final state mass M' = 2.336 to, 
and the same values of other parameters is shown in Fig. [5j For this larger value of the final state mass, the FSI 
contribution is still significant, but it does not dominate anymore. In the real part, FSI and PW contributions 
considerably cancel each other. 

As mentioned above, the form factor F'{Q'^) - which vanishes if the current is conserved - is obtained from 
F{Q^) by the replacement (l23l) in the integrand. The corresponding numerical results for ps = 0.1 and Ps = 0.5 


































11 




FIG. 7: (color online). Transition EM form factor F'{Q^) as a function of for ps = 0.5. Other parameters and 
notations are the same as in Fig. |3] 


are shown in Figs. [5] and [7] respectively. We see that, in comparison with the F{Q‘^) results of Figs. |3][51 the 
value of F'{Q^) is indistinguishable from zero. This very small value is a result of a cancellation between FSI 
and PW contributions. We conclude that in the model considered, the current is conserved. 




FIG. 8: (color online). Transition EM form factor F{Q^) as a function of for ps = 0.1 calculated with the 
’’phenomenological” FSI function (1251) . The notations are the same as in Fig. [S] 


Though the current conservation is natural, the cancellation of FSI and PW contributions is rather delicate 
and provides as a strong test of a calculation. Indeed, both FSI and PW contributions contain the same initial 
bound state BS amplitude. At the same time, FSI contribution contains the scattering state BS amplitude, 
whereas the PW contribution - does not. Their cancellation takes place provided both BS amplitudes - the 
bound and the scattering one - as well as the current operator are consistent with each other, i.e. if they are 
correctly found in the same dynamics. 

To illustrate how a violation of this consistency would affect the current conservation, we replaced the final 
BS amplitude, found for the OBE kernel, by an ” ad-hoc” function 

I 


Tf{ko,z,k) 


(fcp -I- a^){kP‘ + 6^) 


(25) 
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without changing the initial BS amplitude and the current. The transition form factor F{Q^) calculated with 
the function ((25l) for = 1.5, 6^ = 1 is shown in Fig. [8l Apparently it has a typical behavior and nothing 
indicates that it is a wrong result. 

To see that we have displayed in Fig. Othe transition form factor F'{Q^) calculated with the same function 
(USD. We see in this figure that F' is different from zero and of the same order than form factor F{Q^). 
This means that the EM current calculated with the ’’phenomenological” FSI function (1251) is not correct and 
therefore the form factor F{0^) extracted from this current - shown in Fig. [5]- is also incorrect. 



FIG. 9: (color online). Transition EM form factor F'{Q^) as a function of for ps = 0.1 calculated with the ’’phe¬ 
nomenological” FSI function (1251) . The notations are the same as in Fig. [3] 


The study of the numerical stability as a function of the number of Gaussian integration points no, shows 
that the sum F/qj = Fpgj + Fpyp decreases with no- For nc = 64 it is two order of magnitude smaller than each 
Fpgj and F'p-^ taken separately and also than the form factor F. This means that there exist a cancellation 
between F'pgj and Fpyp, and hence the current conservation, with a numerical precision of about 1%. Increasing 
the value of no from 64 to 128 does not improve the result (does not reduce F^^^). Apparently, the precision of 
F/qj is determined by the accuracy of the numeric solutions for the initial and final BS amplitudes. 

The numerical results for F{Q^) and F'{Q^), calculated with na = 64, ps = 0.1 are given in the TableHl For 
= 1.5 the value of F^^^ is by one order of magnitude smaller than F and two orders of magnitude smaller than 
Fpgj and Fpyp). When increases up to = 5, the cancellation becomes worser and almost disappears, 
though F/q( is still a few times smaller than Fpgj and F'p^. This is related to the fact that the BS amplitudes 
was computed in a finite domain of variables k which, at large values of Q^, is not enough to ensure enough 
accurate result for the BS solution. 
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TABLE I: The value F\Q'^) for FSI, PW and tot=FSI+PW vs. in comparison to E(g^); Ps = 0.1. 



E(g') 

F'fsP 

:q^) 

F'pw 


FLtiQ^) 

0.01 

2.948-03 

- 1.571-03 

4.703-02 

-f 

il.732-02 

-4.648-02 

- il.709-02 

5.490-04 

-t i2.304-04 

0.1 

1.391-02 

- il.008-02 

4.541-02 

-f 

H.641-02 

-4.478-02 

- il.633-02 

6.268-04 

-H 18.139-05 

0.5 

1.047-02 

- 7.653-03 

3.904-02 

-f 

il.401-02 

-3.833-02 

- il.383-02 

7.089-04 

-t H.848-04 

1. 

6.640-03 

- 14.841-03 

3.321-02 

-f 

il.181-02 

-3.268-02 

- il.184-02 

5.303-04 

- i3.634-05 

2. 

3.573-03 

- i2.460-03 

2.639-02 

-f 

18.478-03 

-2.537-02 

- i9.150-03 

1.018-03 

- 16.726-04 

3. 

2.362-03 

- il.589-03 

2.234-02 

-f 

i7.062-03 

-2.079-02 

- i7.482-02 

1.552-03 

- 14.204-04 

4. 

1.718-03 

- il.109-03 

1.897-02 

-f 

i6.511-03 

-1.763-02 

- i6.339-02 

1.341-03 

-hil.721-04 

5. 

1.446-03 

- i8.313-04 

1.993-02 

-f 

i6.093-03 

-1.531-02 

- i5.506-02 

4.614-03 

-t i5.873-04 


VII. CONCLUSION 

We have presented the first results of the transition electromagnetic form factor for the electrodisintegration 
of a two-body bound system described by the Bethe-Salpeter equation in Minkowski space. Calculations have 
been performed in a self-consistent way. The initial (bound state) and final (scattering state) BS amplitudes 
were found by solving the equation with the method developed in our previous works |15l | and an OBE kernel. 

We have shown that, provided the bound and scattering state Bethe-Salpeter amplitudes as well as the oper¬ 
ator of EM current are consistent with each other, the electromagnetic current is conserved. If this consistency 
is destroyed, the conservation is violated. This violation has two consequences. 

First, the decomposition of the current in form factors obtains an additional contribution (second term in Eq. 
([3])) which does not satisfy the equality J ■ q = 0. However, the appearance of this term itself does not make 
any influence on observables - it does not contribute in the scattering amplitude - due to conservation of the 
electromagnetic current of the incident electron. 

Second, and most important, is the fact that the non-conservation of the calculated current makes it physically 
meaningless. One cannot extract from a deficient current a reliable transition form factor. It is thus mandatory, 
in practical calculations, like e.g. in the deuteron electrodisintegration, to check the current conservation. If 
the form factor, responsible for non-conservation of the current, turns out to be comparable with the physical 
ones, one can hardly trust the calculated physical form factors too. 

The widely used recipe, consisting in replacing the non-conserved current by the conserved combination 
= Jfj. — Q/iiJ ■ hides the problem but does not solve it. This combination satisfies tautologically the 

current conservation for any J^, not only for the correct one and thus offers no any guarantee to the result. 
With an incorrect current, one cannot find the correct transition form factor, neither from nor from J^. 

The current conservation appears as a numerically subtle phenomenon since it manifests itself as a cancellation 
of large contributions: FSI and PW. To see it unambiguously, the solution of the BS equation should be found 
with high enough precision. 


Appendix A: Kinematics 


As mentioned in Sec. m it is convenient to carry out these calculations in the system of reference where 
p'o — Po, i-e. qo = 0. In the elastic case it coincides with the Breit frame p + p' — 0 and |p| = \p'\, p'q = pq. 

This system exists in the inelastic case M' ^ M too. Indeed, one can easy check that in the reaction 
e + d ^ e' + (np) the momentum transfer q^ = [p — p'Y = (Po ~~ (P~ V'Y i® always negative. In particular, 
its maximal value (reached at the minimal value of s = {M' + meY', with rUe the electron mass) is still negative: 


q^ < — 




M^) 


< 0 . 


M' -I- TO, 
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Therefore one can find a reference frame where this negative value contains only spatial components: = 

{p — p'Y = ~{P~ P'Y so that Pq = po- Moving this frame, without changing the value go = 0) one can also 
make the vectors p and p' collinear, i.e. either parallel or anti-parallel, so that = —q^ = {p — p')^ is equal 
either to (p — p')^ (if p and p' are parallel) or (p-fp')^ (if P and p' are anti-parallel), depending on the value. 
We denote p = | p | and p' =\ p' |. We will precise below when one should take the plus or minus sign. 

In case of an elastic collision M' = M, in the Breit frame with go = 0 (i.e. Po = po) and with the initial and 
final momenta satisfying p + p' = 0 and |p'| = |p|, the scattered system moves in the opposite direction than 
the incoming one and this is the only possibility to get a non-zero momentum transfer. 

In the inelastic case M' > M, still in the reference frame with pg = pg and with the collinear momenta p, p' 
(now with Ip'I ^ IpI but |p'| < |p|), there exists a critical momentum transfer Q^. If the momentum transfer 
is smaller that Q^, it is not enough to change the direction of initial momentum p into the opposite one. 
In this kinematics, the final momentum p' after collision remains parallel to p, though with p' smaller than p. 
However, for the final momentum p' changes its direction relative to p like in the elastic collision. 

That is, when increases, the final momentum p', being first parallel to p, vanishes and appears again in a 
direction opposite to p. When it crosses zero p' = 0 (provided pg = po), the corresponding momentum transfer 
is Ql = p^. We get in this case: 

Po=Po -t yp'^ -b = a/p^ -I- M2 
^ M' = vp^Tm^. 

From last equality we find the critical value = p^ for which p' = 0 

Ql = - M^. (AI) 


In the frame where we perform the calculations (go = 0), the following kinematical relations holds: 

(A2) 
(A3) 
(A4) 
(A5) 
(A6) 


(A7) 

and denote hereafter (abusively) : k =\k\. 

From the requirement a/ M'^ + p'^ = A^M^'^rp^ we find relation between p and p': 

p' = 

p = \Jm'^ - M2 -Kp'2 

and so: 






{p'-pf = {p'-apf 
k ■ p' = cr kp'z 
{p — kY = p^ — ^zpk + k"^ 

(p' — fc)^ = p'"^ — 2azp'k + k^ 

= p- up 

where we have introduced the ’’sign” variable u depending on and Q\\ 


<y{Q\Ql) = 


+ I Q2 < q2 
-I ij Q2 > Ql 
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P 

Po =Pq 


Ql + Q^ 

2 v^ 


Q 2 )((M' + M) 2 +Q 2 ) 


2^ 


(A 8 ) 

(A9) 


Appendix B: Calculating the FSI contribution to the transition form factor 

The contribution of FSI to the transition form factor is given in Eq. (jl 6 p as a sum of three terms F 3 2 ,i- In 
their turn, -^ 3 ^ 2,1 are obtained by integrating over dk^d^k the three functions /a, 2,1 defined in (1I3II51). We detail 
in what follows the calculation of these three contributions. 


1. F3 contribution 


Let us first consider the contribution 


F^{Q^) = 


(27r)3 


roo n+l 

J J ' 


^ + 00 


dz dko f3iko,z,k) 


with /a given by Eq. 0 . 

As a function of /cq) fs contains six poles in fco-variable at the points 

ko = FSj: 
ko = -Sj: 
ko = Po+ep-_fe 
ko = Po-£j,_u 
ko = Po+£p_k 
ko = Po — £p-k 

and eventually other singularities resulting from G(ko, z, k). We subtract and add to /a a function /13 depending 
on the same variables and having poles only in the variable kg. That is 

/s = /a — ^3 + ^3 = /a + ho (Bl) 


with 


h3{ko,z,k) 


gi{z,k) ^ g2{z,k) ^ gojz^k) ^ 9^{z,k) ^ 9b{z,k) 

fco-£fc ko+e^ ko-po-£p_k fco-Po + £p._fe ko-p'o-£p_i: 


9&{z,k) 

ko-p'o + £p-k 


(B2) 


The coefficients gt, independent on ko, are determined by imposing that the difference 


fs = fa — ho 


(B3) 


is regular in fcg- For instance 

lim (fco - £k)fo = 0 lim {ko - £k)fo = gi{z, k), 

ko — ko—>Sk 


and similarly for other poles. 
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This gives 


9i(z,k) = 
92 {z,k) = 
93{z,k) = 

94{z,k) = 
95{z,k) = 

96{z,k) = 

and the function /a obtains the form: 


G{+£k,z,k) 

+ 2e,[ipo - Sk? - el_^] [{po - Sk)^ - 

_ G{-ek,z,k) _ 

2efe[bo + £k)^ - el_k] [{Po + £k)^ - e2,_^] 

G{po + ep-k,z,k) 

+ 2ep-k[ipo + £r-kr-£l] 

Gi^Po £p—k-i Zjk^ 

2ep-k[{po - £p-kf - e|] (^p-fc - £l'-k) 

_ G{po +epi-k,z,k) _ 

2ep'-fc[(po + £p'-kY ~ £fe] {£p-k ~ ^p'-k) 

G{po- ep'-k,z,k) 

^ 2ep,.u[{pc-£p'-kY-£l] {£U-^l-k) 


(B4) 


h{ko,z,k) 


_ G{ko,z,k) _ 

[ko - £k] [ko + £k] [ko - Po - £p-k] [ko - Po + £p-k] [ko 

9i{z,k) _ g2{z,k) _ 93{z,k) _ 94{z,k) 

kg — £k ko + Sk ko — Po — £p-k ko — po F £p-k 


~ Po ~ £p'-k] [ko — Po F £p'-k] 

_ 95{z,k) _ 9o{z,k) 

ko Po £p'—k ko Po “t“ £p 


(B5) 


' —k 


The PV integral over fcg of the remaining integrand /13 in 


vanish in the full integration domain: 


PV 


ho dko = 0. 


(B 6 ) 


However, in the numerical solution we restrict the integration domain to a finite interval ko G [—L,+L]. The 
integral (IB 6 I) is no longer zero and finite volume correction must be taken into account. The integral over the 
finite domain [—L, FL] of the function ho given in (IB2p is analytic and reads: 


/ -r-G 

dko ho{ko,k,z) = gi{k,z)log 
F go{k,z)\og 
F g 5 {k,z)\og 


L-£k 


L F £k 

T Po £p—k 


L F Po F £p—k 
L — Po — £p'-k 


L F Po F £p'—k 


Fg2{k,z)\og 
■f 54 (fc, 2 ) log 
+ S6(fc,^)log 


L F £k 


L-£k 
L ~ Po F £p-k 


L F Po ~ £p—k 
L — Po F £p'—k 


L F Po ~ £p'—k 


For the raisons that will become clear latter, we will include the above finite volume contributions 


= 




(27r 


dz fojv{z,k) 


(B7) 


(B 8 ) 


in the F 2 contribution, to be discussed in the next section. 

The Fo{Q^) will thus be given by the following three-dimensional integral: 


FoiQ^) = 


(27r)^ 


k^dk / dz dko fo{ko, z, k) 


(B9) 
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with /^(koyZ, k) defined in (IB5I) . 

The integrand in (IB9I) is by construction a smooth function in fcg. Concerning the integration over the 
variables k and z, a further inspection is needed. For M' > 2m the value po — eg — e^,_g = p'^ — 
in the denominator of (IB5I) can vanish. Indeed, a state with mass M' > 2m can decay in two particles with 
masses m, and in the moving frame this implies Pq = eg + The gi and ge terms in Eq. (IB5I) are therefore 

singular in variable k. These two singularities cancel each other since an expansion of /a in the variable eg near 
^k~P^~ ^p'-k shows that the result is proportional to ^ C>((e^,_g — Po)°)' ^s a consequence, the integrand 
/a has no additional singularities in fc, z and can be safely integrated numerically by standard methods. 


2 . F2 contribution 


Let us now calculate the contribution 


F2{Q^) 


i 

(27r)S 



f2{ko,z,k) 


(BIO) 


with /2 given by Eq. (HI- 

The integration over fco"Variable can be performed analytically by means of the delta-functions. The result - 
denoted /2 - is expressed in terms of functions gi defined in (IB4I) and reads: 

/ +00 

dko f2{ko,z,k) = -iir 

-00 


[giiz, k) - g2iz, k) + gsiz, k) - g^iz, k) - 55 ( 2 :, k) - p 6 (^, fc)} (Bll) 


As one can see, f 2 {z,k) has a similar structure and depends on the same variables than the finite volume 
corrections f^jv{z,k) described above in (IB7I) . It is thus natural to include both contributions in the same 
integrand 


6 

Mz, k) = hiz, k) + hjv{z, k)=Y^ Ci{k)g,{k) 

2=1 


by introducing the coefficients 


ci(A;) 

= —iiT 

+ 

log 

C2{k) 

= Fiir 

+ 

log 

csik) 

= —in 

+ 

log 

C4{k) 

= +in 

+ 

log 

Cb{k) 

= +in 

+ 

log 

ce{k) 

= +in 

+ 

log 


L-Sk 

L+Sk 

L-\-ek 

L-€k 

L Po Sp — k 

L-\-po-\-£p — k 

L—po-\-ep-k 
L-\-pQ — Sp-k 
L-po-€^f_k 
L+pq+Sj^/ _k 

L-po+£p/-fc 

L+PQ — £pr _k 


(B12) 


The F 2 contribution is then given by the two-dimensional integral 


F2{Q^) 


r+1 


(27r)3 


k^dk / dz f 2 {z,k) 


'-1 


(B13) 


The integral (IBlfij) over the z variable requires some care since both gi and g^ can have pole singularities in z. 
In contrast to the function /a, Eq. (IB5I) . these singularities in /2 do not cancel each other. However, they can 
be integrated analytically over z, so the pole singularities turn into the log-ones. 
















2 F 2 contribution 


18 


Let us first consider the gi term: 


9i(z,k) 


G(+Sk, z, k) 

2e,[{p,-euY-el_,] [{p, - 


The denominator vanishes if 

(Po - - Ep'-fc = 0 2poek-M'^ = 2p'■k = 2ap'kz 

with CT = ±1 is the sign function defined in (IA7|) . A singularity in the 0 -variable would exist a.t z = zq given by 


zoik) = a 


2po£k - M 
2p'k 


/2 


provided | zq |< 1, that is for k in the interval k- < k < k+ with 

1 


k^ = - 


p'TPo\ll-{ ^ 


(B14) 


(B15) 


Notice, in particular, that this singularity exists only in the inelastic case since 2m < M'. It is a moving 
singularity, depending on the value of the second argument k, as well as on the momentum transfer and the 
parameter Q^. 


To properly account for this singularity we split the ^-integration interval in three domains 

[0, -boo] = [0, kJ\ U [k-, fc+] U [fc+, -boo] 
as well as the corresponding integral (IB13I1 


F 2 {Q^) = I1+I2+I3 


with 

■ ^- 1-00 ^- 1-1 

^3 = / dz f 2 {z,k) 

[2'^) Jk+ J-1 

• The integrals over [0, fc_] and [fc+, -boo] have a smooth integrand in both variables and the contributions 
Ii and I 3 can be computed by standard methods. 

• The integral over [fc_, fc+] has a singularity on at z = zq. The integrand is regularized by using the usual 
subtraction procedure: 


9 i{z,k) 


9iiz,k) 


9'i{k) 

Z- Zo 


9[{k) 
Z- Zo 


(B16) 


g[{k) = lim (z - zo) 5 i(z, k) = Res [ 5 fi(z, k)] 

Z—^Zn ^ 


with 
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To compute this quantity we take z from 

Cp'-k = + p'^ — 2 ap'kz = + Pq — — 2 ap'kz 

and together with (IB14I) we have 


4 + Po - 9-fe - 2po£fe bo - Skf - cl,-k 
z — Zq = a -- = a- 


el+pl-M'^-el, 


2 p'k 


2 p'k 


and so 


G(+£fe, k) 


(z - z^)g,{z, k) = a— _ , .2 _ .2 


We get in this way the residue 


Ap'kek[(po - euf - el_^] 

(lA - ^ f G{ek.z,k) 

^ 4p'fc£fe I 


bo 4 “ ^p—k) bo ^k ^p—k} J jg— 2 q 


2 p'k 


The ge term 


96{z,k) = 


G{po - ep>-k,z,k) 


2Spf — kipc ^p' — k 4” £^fc)(P0 ^p' — k ^k) i^p—k 4“ Sp' — k)(^Cp—k ^p' — k) 

has the same singularity at z = zq than g 2 and has been treated in the same way by subtraction 


96iz,k) = 


96{z,k) - 


g'eik) 


z - Zo 


+ 


g'eik) 

z - Zo 


By a similar calculation we find 


g'eik) = -g'lik) 


Finally, once regularized the singular terms g 2 and ge, the I 2 contribution is given by two integrals, 
spending to the two terms in the subtraction (IB16|) 


I2 = H2+I2 
i 


/’+^ 

k^dk J dz f 2 {z,k) 


T" — 
-'4 — 


( 27 r) 


, k dk c[{k)g[{k) log 

Iki 


1 - Zo 


1 4- Zg 


where 


f' f 

12 = / 2 -cbfc)- 


z - zo 


(B17) 


corre- 

(BI 8 ) 

(B19) 

(B20) 

(B21) 


is a regular integrand, g[ the residue (IB17I) and Ci{k) = ci(fc) — CQ^k) with Ci given in (IM 
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3 . Fi contribution 


The Fi contribution is given by 

Fi(Q^) = ^^^^3 J dko dz k'^dk fi{ko,z,k) (B22) 

According to definition 03) of /i, this function contains three contributions, each of them contains the product 
of two delta-functions. It turns out that the non-zero contribution results from the second term only. 

We substitute m in (|B22I) . integrate this term over ko analytically and find: 


fiiz,k) = dkofi = - 


Gi£k.z,k) 


4ere-/ t —e 

k p'—k ^ — h 


S{po - eg - ep_f:) 


p' — k 


p—k' 


(B23) 


to be integrated then over z and k: 


FiiQ^) = I dz k^dk /i(z, k) (B24) 

In general, the denominator (eg ~ e?. -) in (IB23P can vanish. However, one can check that if the argument 
of the delta function (pn — st — s t) is zero - providing a non-zero contribution - the denominator (e - - —e?, 

' K p —K' ± cj \ pf_]^ p—k' 

is not zero. Therefore there is no singularity from this denominator. 

Concerning the first and the third terms in (jI5p . the arguments of the (5-functions in them vs. ko can also 
cross zero. However, in the first term they cannot be zero simultaneously. In the third term, the arguments 
of the (5-functions can be zero simultaneously. That is, after integration over ko, we obtain the delta-function 
^ S{E^ — A?. -) which could contribute. However, a more careful analysis shows that its contribution is in 

p—k p' — k 

fact zero. For this aim we represent the delta-function as 


6 {x) 


e 

7r(x^ -I- e^) 


and integrate over both dz and dk. Taking after that the limit e —>■ 0, we find a zero result. Care must be 
taken however to do not take this limit too early, i.e., before integration, since we will get in this way a wrong 
non-zero contribution. 


After integrating analytically over dz in (IB24I) by means of the delta-function we obtain: 

/: 


Mk) = / A* = - 1.„|) 


p' — k' 


(B25) 


2 = 20 


where the value of zo is given by (IBI4I) . 

Finally, the integration over k is reduced, due to the theta-function 0(1 — l^ol), to the interval k G [fc_,fc+] 
with kzfi given (IB15I) . That is: 


FiiQ^) 



k^dkfiik) 


(B26) 


As a test, we carry out an independent calculation for Fj = Fg = 1, using Feynman parametrization for the 4D 
integrals. In this way, we find the imaginary part which coincides with the contribution (jB26l) . 
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4. About the function G{ko,z,k) 


An important quantity used in our formalism, which contains all the information about the initial and final 
state vertex amplitudes, is the function G(/co, z, k), defined in Eq. dill) 


G{kQ,z,k) = Ti - k,pj Tf 


Some useful relation concerning this function are specified in what follows. 

The initial state amplitude ri(A:oj k) — we denote \k\ = k — is computed in the reference frame where p = 0. 
In an arbitrary frame Ti(ko, k) is written as ri(fco, k) where 

r k-p 
kn = — 


M 


k = 


{k-pY 

M2 




The arguments of T^ (| — fc,p) are obtained from these relations by the shift fc —>■ ^ — k. The same happens for 
the final state vertex amplitude T/ — k,p''^. The function G(fcoj -Z) k) should be therefore understood as: 


G{ko,z, k) = ^r,(fco, fc)r/(fc'o, k') 


Po 


(B27) 


with, after performing the shift k ^ ^ — k, the arguments are given by 


kn = 


k = 


M k-p 

Y~ Jf 


{k-pY 




Written in more detail: 


kn 

k 

k'n 

k' 


M2 

M knpn — kpz 

T ~ M _ 

_ l {koPo-kpzY 2 r, 
- V M2 ^-^o+k 

M' knPo — k-p' 

~ Y~ M' _ 


M'^ 


fc2 


(B28) 


We remind that the sign in the scalar product k-p' is given by 
by Eqs. dA8l) and dA9|) . They all depend on the value of 


and the values of p, p' and po are defined 


The initial (bound state) solution T^ is normalized so that the elastic EM form factor at = 0 is 1. There 
is no any uncertainty in the normalization of the scattering state solution T f de termined by the inhomogeneous 
BS equation. One should also take into account that the solution found in [15l | was the partial wave amplitude 
Fq related to the full amplitude by 


Tf = 167r^F,Pi(cos6l) 


;=o 


and that for the S-wave, the function T/ in (IB27I) is related to our solution Fq obtained in [T^ by 


T/ = IGttFo. 


(B29) 
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Appendix C: Calculating the PW contribution to the transition form factor 


We carry out here the part of integration over and d^k in the form factor (HI that can be done 

analytically. Integrating first over dfljs^ in the frame p' = 0, pg = M' = 2^Jrn? + pi we get: 




df^p^ 

47r 


1 


47rp2 


<5 Uo- 


M' 


d(|fc| - IpII) 


In an arbitrary frame ko and | k \ are rewritten as 

k-p' 


M' 


\k\ = ^kl - k^ 


jk-p') 

M'2 


-P 


Explicitly, in the frame where pg = po'. 


k-p' kopo — az\k\p' 

'M' ~ W' 




where the sign a is defined in (IA7I) . 

After these transformations, the form factor Eq. (ITS)) obtains the form: 

(po - fcg) r» (I - k,p) 

Po [(p - fc)2 -rn? + ie] 

([koPo - (^zkp' M'\ f I (fcoPo - akp'y 

" H— M' -— IP — 

Eq. (ICll) contains two-delta functions and integration over three variables fcoj k and z. By means of the delta- 
functions we integrate over the variables fco ^md z. If argi and arg 2 denote the arguments of the first and second 
delta-functions, the conditions that the arguments of both delta’s equal to zero give the system of equations 
argi = 0, arg 2 = 0 which we solve relative to fco and z. We find that ko = Sk = \/m? -|- fc^ and z = zo with the 
value of Zo given by (IB14I) . Remind that | zo |< 1 if ^ G [k-, fc+] with k^ given by (IB15I) . 

Calculating an integral containing the delta-function, one should divide the result by derivative, over the 
integration variable, of the argument of the delta-function. Since the integral (ED contains two delta-functions, 
we have to divide the result by the product di d 2 of derivatives di = arg '^, d 2 = arg '2 of the arguments of both 
delta-functions. Each derivative depends on the order of calculation (first over fco, then over z or in the opposite 
order), though the final results, which is determined by their product, is the same. That is: 

77 / / II kskP 

di d2 = arg^^ko 0-^92,z = ar9i,z arg2,ko = 


fc^ + fc2 - ps 


dkok^dkdz 

2 ^! 


(Cl) 



provided argi = arg 2 = 0. 

Einally, the integral over z is reduced to: 



... 5{z — Zo)dzdk 


...9{l-\zo\)dk = 



...dk 
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Thus, after integration over ko and z we obtain the integral (1191) over k in the limits determined by the condition 
\zo\ < 1 and given by Eq. (jB15|) . 
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